Dissipative Quantum Ising model in a cold atomic spin-boson mixture 
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Using cold bosonic atoms with two (hyperfine) ground states, we introduce a spin-boson mixture 
which allows to implement the quantum Ising model in a tunable dissipative environment. The first 
specie lies in a deep optical lattice with tightly confining wells and forms a spin array; spin-up/down 
corresponds to occupation by one/no atom at each site. The second specie forms a superfluid 
reservoir. Different species are coupled coherently via laser transitions and collisions. Whereas the 
laser coupling mimics a transverse field for the spins, the coupling to the reservoir sound modes 
induces a ferromagnetic (Ising) coupling as well as dissipation. This gives rise to an order-disorder 
quantum phase transition where the effect of dissipation can be studied in a controllable manner. 

PACS numbers: 03.75.Mn, 64.70.Tg, 71.27.-fa 



Spin-boson models are essential ingredients in quan- 
tum optics [l| , nuclearphysics [4I , quantum chaos |3| , and 
quantum dissipation [J, |5|- In particular, an ensemble of 
identical two-level systems, each coupled to common ra- 
diation field modes, the Dicke model |6|, was introduced 
initially to describe the superradiant emission (a sudden 
increase in the rate of coherent spontaneous emission of 
an ensemble of atoms). This model exhibits an inter- 
esting order-disorder transition for the two-level systems 
both in the limit of only one or several boson modes [7| . 
In this Letter, we theoretically envision a different "spin- 
boson" mixture, i.e., a spin array coupled to a large col- 
lection of harmonic oscillators, realized using cold atomic 
bosons, which allows us to explore how the properties of 
the celebrated quantum Ising model and of the emergent 
quantum phase transition for the spins [8| are modified in 
a tunable bosonic environment. We emphasize that our 
setup embodies the first tunable realization of the quan- 
tum Ising model in a dissipative bath, and that several 
critical exponents can be measured using standard imag- 
ing techniques. Generally, mixtures of different species of 
cold atoms open a fascinating field of many-body physics 
where exotic spin systems and quantum phase transitions 
can be engineered and probed [9|, [lOJ, |ll | . 



More precisely, we consider cold bosonic atoms with 
two (hyperfine) ground states a and b, trapped by dif- 
ferent state-selective external potentials [l2[. Whereas 
atoms in state a form a Bose-Einstein Condensate (BEC) 
of dimension d — (1,2,3), held in a shallow potential 
14 (x), the h atoms are trapped in a deep d-dimensional 
optical lattice with well-separated lattice sites, so that 
hopping between different sites can be neglected (see 
Fig. [1]). We consider the coUisional blockade regime of 
large onsite interaction Ubb, where only states of occupa- 
tion number nf, = 0, 1 per lattice site contribute to the 
dynamics. With all higher occupied states being adi- 
abatically eliminated, one obtains a d-dimensional ar- 
ray of two-level systems or quantum 1/2-pseudospins. 
Spin-up/down corresponds to an occupation number of 
one/zero. The spins are coupled to the low-energy 



FIG. 1: Setup comprising a BEC and two-level systems. 



phonon excitations (superfluid sound modes) of the BEC 
in two ways: via collisions between a and b atoms and 
via laser mediated Raman transitions between the two 
hyperfine species. The dynamics of a single two-level 
system embedded in a BEC has been addressed in Ref. 
U and the emergent quantum phase transition can be 

understood in terms of the spin-boson model [4, H , which 

• \\ I — li 

has abundant practical applications in physics [141 llSj . 

The coUisional interaction of the atoms is described 
by a contact pseudopotential with coupling strength 
ffa/3 = 'iTraaish^/m for a, /3 G {a, 6}, where Uajs is the cor- 
responding s-wave scattering length and m the atomic 
mass. As shown below, the cross-coupling gat between 
the spin array of b atoms and the superfluid reservoir of 
a atoms induces a ferromagnetic (short-range) Ising ex- 
change interaction between spins. On the other hand, at 
each site, the Raman coupling acts as a transverse field 
A, favoring paramagnetism. In addition, the coupling 
of the spins to the long-wavelength BEC-phonons intro- 
duces dissipation to the dynamics of the spin system. In- 
tegrating out the BEC (phonon) degrees of freedom, we 
demonstrate that this will result in a dissipative quan- 
tum Ising model (DQIM) for the spins [3, 111, [3, that 
exhibits an unusual dissipative order-disorder phase tran- 
sition. One can determine the phase of the spin system 
(ferro- or paramagnetic) by measuring the mean occupa- 
tion number {ni,) of the lattice sites [l9|. In addition, 
measuring (rif,) for different Raman laser intensities and 
detunings allows to extract several critical exponents. 

Model.- The a atoms form a BEC (which can have di- 
mension d = 1,2,3), and the Bose-ficld operator can be 



split into magnitude and phase V'a(x) = ^/5a(x)e~"^''''. 
We consider the low temperature limit where the su- 
perfluid density is equal to the equilibrium atomic den- 
sity Pa- Introducing the density fluctuation operator 
n(x) — Pa{x.) — Pa, which is canonically conjugate to 
the superfluid phase c^(x), the dynamics of the BEC 
in the long-wavelength approximation is then described 
through a quantum-hydrodynamic Hamiltonian, 



Ha 



dx 



2m 



|V0(x)|^ 



raVg 



ifU) 



(1) 



where Vs is the sound velocity and m the atomic mass. 
This description can be applied to all dimensions d = 
(1,2,3). In d = (2,3), one obtains this Hamiltonian re- 
sorting to the Gross-Pitaevskii approach where the in- 
teraction coefficient gaa enters the Hamiltonian via the 
sound velocity v1 — PaQaa/'m. In the one-dimensional 
(ID) case of a tight transverse confinement with fre- 
quency fiw±^ 3> (ksT^p) (with /i being the chemical 
potential), which is particularly interesting as it leads 
to ohmic dissipation, the proper derivation follows the 
Haldane-Luttinger approach [2Q] . In fact, one may iden- 
tify K = hTrpa/{mVs) as the Luttinger parameter, where 
Pa = Pal^ is the ID density that depends on the trans- 
verse confinement length l± = {h/mu!±)^''^ . The corre- 
sponding ID interaction coefficients are related to the 3D 
scattering lengths by (jafj — 2huj±aai3 if Uap <^ l± |21l] . 

The Hamiltonian Ha can be diagonalized using phonon 
operators via 0(x) = z^q I 2l^^l'^'e*''"(7q-7iq) and 

n(x) = Eq I 2^l'/'e^''''(7q + 7lq), whcrc V^L'^is 
the volume of the (BEC) system. Since the low-energy 
excitations of a BEC are phononlike only for wavelengths 
larger than the healing length ^ = h/y/2mpagaa, the sum 
is cut off at a frequency uic = Vs/S, = \/2gaaPal^- One 
gets the well-known commutation relations, 

[0(x), n(x')] = f E ^'^-''' = » "^'";f'^^^ (2) 

q 

where each spatial component of the wavevector runs be- 
tween L~^ < \qi\ < ^~^, and ?7(|x— x'|/^) is a function 
which decays very rapidly for separations larger than 
the healing length. For instance in d = 1 one finds, 
r]{xjk) ~ sm{xjk)/xjk with Xjk = |xj - xa;|/^. Note that 
for C — > 0, one checks that [4>{x.), II(x')] -^ i6{x — x')). 

The Hamiltonian (Hb+Hab) that describes the dynam- 
ics of the b atoms as well as the coupling between the two 
species (colfisions and the Raman coupling), reads [l3|: 



Hb + Hab = 2_^ 



-hro+ / dx\lJjbAi-x)\'^gabPa{^) 



Hb^ 



iyitel [V;a(x)V'b,,(x)&|-|- 



h.c. 



Ubbb]b]bA 



Here, b] creates a particle in the localized Wannier 
ground state V'fa.i at lattice site i which has a spatial ex- 
tent lb- The operator V'a(x) annihilates an atom (in state 
a) of the superfluid reservoir at position x. Transitions 
between the two hyperfine states are induced by Raman 
lasers with Rabi frequency f2 and detuning rg. Atoms 
are transferred from the reservoir to the lowest vibra- 
tional band of the lattice and vice versa. In coupling the 
b atoms to a coherent matter- wave field ■!/'a(x), we have 
assumed that the number of a atoms inside a lattice well 
is large, Ua = Palf ^ ^- We also include the colfisional 
interaction between a and b atoms. The repulsive onsite 
interaction Ubb ~ 9bb/lb > of & atoms can be varied 
using a magnetic or optical Feshbach resonance [l3| . We 
assume Ubb/h ^ (^jJ'o) in order to neglect doubly and 
higher occupied states (colfisional blockade limit). 

Using Pauli matrix notations to describe the spin ar- 
ray, we replace the occupation number operator &J&i by 
(1-t-crf )/2 and bj —i- af . Furthermore, in the spin-phonon 
coupling term, we can replace the phonon field operators 
by their values at the center of the lattice sites x^ , because 
the phonon wavelength is much larger than the extent of 
the Wannier ground state lb- Neglecting constant terms 
we obtain the Hamiltonian, 



Hh + H„ 



E(-? + ^n(x.)).r 




(4) 



1/2 

where the sum runs over all lattice sites and A ~ Ua f2 
is an effective Rabi frequency which grows with the 
number of a atoms inside one lattice well. This form 
of the Rabi coupling requires that the cutoff frequency 
^c ~ Padaa ^ ^A, siucc wc Only consider "single- 
particle" hopping events. We have also introduced the 
effective detuning hr — hr^ — pagab + {^^Y /{'2'Ubb) which 
includes a colfisional mean field shift as well as a shift due 
to the virtual admixture of doubly occupied states. 

Phonon integration for the mixture. -To find the dy- 
namics of the spin array, we shall integrate out the 
phonon modes along the lines of a single two-level system 
embedded in a BEC [13]. In the case of a spin array cou- 
pled to bosons in the superfiuid phase, we show that the 
integration of phonons can be still carried out exactly, 
which is in contrast to the case of a spin array coupled 
to conduction electrons (Kondo lattice) 18]. 

Firstly, we apply the unitary transformation H' = 



WHU with U = exp -il]^ cr|^(xj)/2 to the total 
Hamiltonian H = {Ha + Hab + Hb) in order to eliini- 
nate the operator exjp{—i(f)(xi)) from the Raman term. 
As a result, the spin-bath coupling takes the exact form. 



(3) 



E 



-1 



9ab 
9aa 



hA 

2paga. 



ffaan(Xj) 



(5) 



where the cr^-term can be neglected as long as one stays 
away from the point gab = gaa since A <C Wc- As we 
show below, we can only reach the interesting regime 
of the quantum phase transition if \gab\ ^ gaa, since the 
ferromagnetic interaction dominates otherwise. Thus, we 
assume that gat is sufficiently small and keep only the 
coupling term of the z-component of the spin to the bath. 
The unitary transformation also produces a spin-spin 
exchange term, and the total Hamiltonian becomes. 






ha^ 



E 

3,k 



PagaaVjk 



1 ^9ab , . , 

9a ' 



E ^"^iT'q 



7q: 



(6) 



where we have written ?7(|xj — Xfe|/^) = rjjk and o^q = 
iis|q|. For realistic conditions, the separation between 
neighboring lattice sites db is of comparable size as the 
heahng length {(, ^ db > h) [l3|, and since the function 
rjjk defined in Eq. ^ decreases very rapidly for separa- 
tions larger than S,, below we can safely restrict ourselves 
to nearest-neighbor spin interactions. We have defined: 



K 



9abPa 



mvt 



1 



mv^q 



2Vpah 



1/2 



(7) 



Now, the integration of phonon modes can be carried out 
exactly (using coherent state functional integrals). This 
produces an extra term in the action of the spin array, 



SS' 



^q 



E 



drdr'-^ 



icOqiln 



xe^""(^-^'Vq(r)<(r'), (8) 



where /3 = (fc^T)-!, a^ir) = E, ^|(^)e^''"^ ^n = 
2im/{h[3) are Matsubara frequencies, and Aq/wq does 
not depend on q. In general, the coupling to the bo- 
son bath provides two distinct contributions, which can 
be identified through, i7q(r)i7q(T') — ^(o'q(T)crq(T) -f- 
^q(T')<(r')) - W^ir) - aq(r'))(a*(r) - a*(r')) [5|. 

The term which is local in time simply renormalizes 
the spin-spin Ising interaction. Adding up all spin in- 
teraction terms, we find that the overall interaction is 
ferromagnetic and is limited to nearest neighbors (as the 
healing length ^ « db): J2{j.k) [•^/'^'^ '^|('^)^fc('^)] '^ith 
J = —gab/{^'^9aa£,) < 0. In coutrast to itinerant elec- 
trons which produce an RKKY interaction between spins 
which can be either ferromagnetic or antiferromagnetic, 
dependently on the electron density, here phonons result 
in a ferromagnetic Ising exchange [23|. We have checked 
this result by setting A = and by applying a different 
unitary transformation U' = exp[— i^, ^^^^ '^j 4'i'^j)] to 
the Hamiltonian H in accordance with Refs. 122, l23l. 



Similar to the case of a single two-level system coupled 
to a bosonic bath [5|, |24| , the second (dissipative) contri- 
bution stems from the pole at ujq — — iJ7„ (or ujq = ifln). 
Finally, the spin dynamics is dictated by the action [241] . 



2 ^^ ' 2 ^^ ' 



-gE^ dTdT'a{T-T')a]{T)a]{r' 



(9) 



a(T — t') = i Jp°° (icjJ(a>)e^'^l'^^'^ ' involves the spectral 
density J(w) = fm Eq |Aqp5((j-' — Wq) = 2TThauj'^ (iJ7„ -^ 
Lu + iO^) and the dimensionless dissipative parameter 



Do{d)mv'^ 
4:h{27ryvfpa 



9abPa 



mvi 



(10) 



where Do(l) = 2, Dq{2) = 2tt, DoiS) = in and in the ID 
case one must use the variables pa — > Pa and gaf3 — > ga(3- 
Hereafter, we focus on the limit r ^ A, where one recog- 
nizes the d-dimensional DQIM of Ref. [16| . The Raman 
coupling A favors a paramagnetic state where (erf) = 
which competes with the ferromagnetic exchange J. The 
effect of dissipation will be more prominent in the ID 
case, which corresponds to ohmic dissipation (d = 2, 3 
rather corresponds to super-ohmic dissipation) [J, |5| . 

Theoretical Predictions.- For gab = 0, the system will 
be in a paramagnetic regime (since J = 0) , in accordance 
with Eq. (4), whereas for large \gab\ '^ 9aa > ^^1 Pa, 
the system will be in a ferromagnetic phase. One can 
determine the phase (ferromagnetic or paramagnetic) of 
the spin system by measuring the mean occupation num- 
ber [ub) — ((erf) + l)/2 19]. The dissipation produces a 
damped Rabi oscillation in the paramagnetic phase [l3|. 
Below, we will assume \gab\ ^ 9aa in order to study 
the emergent "dissipative" ferro-paramagnetic transition; 
this results in < a < 1 (in ID: a « l/(4i^) < 1/4). 

The nondissipative quantum Ising model is known to 
exhibit a second order phase transition around | J| r^ ftA, 
separating the paramagnetic and ferromagnetic phase, 
and can be mapped exactly onto the classical Ising model 
in dcs = d -\- z = d -\- 1 dimensions; z is the dynamic 
critical exponent that describes the relative dimensions 
of (imaginary) time and space. The case z — 1 implies 
that time just acts as another spatial dimension. The 
classical model itself can be described by a ^"'-theory [8| 
with well-known critical exponents [25| . 

In particular, the quantum Ising chain lies in the uni- 
versality class of the 2D classical model 26[ . With dissi- 
pation (a > 0) , the critical behaviour is most profoundly 
changed in the ohmic case {d = 1). The bath not only 
renormalizes the transverse field. A' = A (A/wc)" ~" 
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FIG. 2: ID phase diagram as a function of a and J. 

[J, |27|, such that the transition occurs at a smaller 
\J\ ~ fiA' (see Fig. [5]), but it also changes the univer- 
sality class of the quantum phase transition. The critical 
exponents of this second order phase transition are inde- 
pendent of the value of a [8|, ll6|, [iTf and they have been 
derived through a d issip ative ^^-theory in d = (2 — e) 
dimension (e = 1) [l7|, |l8| and through Monte Carlo 
simulations [l6[. Below, we summarize the main results. 
The dynamic critical exponent is equal to z = 2 in the 
ID case, whereas z = 1 for d = (2,3). Thus, the ID 
and 2D dissipative quantum Ising models should behave 
similarly, since they have the same effective dimension 
dcff = d+ z Ki 3. Experimentally, one can measure [8l.l25| 
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(11) 



dro 



where (erf )c refers to the value of (trf ) — 2(6- 6^) — 1 at the 
quantum phase transition in the presence of a finite de- 
tuning ro and Ac refers to the value of the Rabi frequency 
at the transition. The exponents can be extracted us- 
ing standard state-selective measurements of the h atom 
mean occupation number, averaged over an homogeneous 
part of the lattice. By measuring the three exponents in 
Eq. (11), one may confirm that the phase transition is of 
second-order type where they are related by /3((5— 1) =7. 
In the case of the ID dissipative system, the value of 
the dynamic critical exponent can be directly obtained 
from the equality z = {b — l)/2 ~ 2. In addition, from 
scaling laws, one finds that 7 = zz/ « 1.276 [H, O El 
and /3 — v/2 « 0.319, which implies that the correlation 
length exponent v « 0.638 can be inferred as well. Note 
that (erf) goes continuously to zero at a second-order 
phase transition which is in striking contrast to the case 
of a single two-level system embedded in a ID BEC 13 1 



where a Kosterlitz-Thouless type transition takes place, 
and where (rif,) shows a finite jump as the detuning is 
changed from a negative to a positive value [l5| . Table H] 
summarizes the critical exponents for various systems. 
In short, we have shown that a spin-boson mixture of 
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Dissipative d = 1 


0.319 
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Non-dissipative d = 1 


1/8 


7/4 


15 
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(Non)-dissipative d = 2 


0.325 


1.241 


4.82 


0.630 


1 



TABLE I: Critical exponents of quantum Ising models. 



cold atoms can be used to engineer the quantum Ising 
model in a dissipative bath. System parameters can be 
tuned to reach the quantum critical regime and critical 
exponents can be measured. Finally, the spin dynamics 
of the DQIM remains an open and challenging question. 
Discussions with P. ZoUer are acknowledged. This 
work was supported by NSF through the Yale Center 
for Quantum Information Physics. 
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